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Abstract. We study a particle model for a simple system of partial differen- 
tial equations describing, in dimension d > 2, a two component mixture where 
C^ ■ light particles move in a medium of absorbing, fixed obstacles; the system 

consists in a transport and a reaction equation coupled through pure absorp- 
Wj ' tion collision terms. We consider a particle system where the obstacles, of 

radius e, become inactive at a rate related to the number of light particles 

travelling in their range of influence at a given time and the light particles are 

' (^ ' instantaneously absorbed at the first time they meet the physical boundary 

^ ' of an obstacle; elements belonging to the same species do not interact among 

Q ^ themselves. We prove the convergence (a.s. w.r.t. the product measure asso- 

O ■ ciated to the initial datum for the light particle component) of the densities 

describing the particle system to the solution of the system of partial differ- 
ential equations in the asymptotics a'^n~'^ — > and a'^e'^ — >■ 0, for re g (0, n) 
►-^ , and C € (0, 2 — 2d)' "here a^ is the effective range of the obstacles and n is 

VO ■ 1- Introduction 



the total number of light particles. 



In this paper we propose a microscopic model for a system of partial differential 
equations consisting in a transport equation, having pure loss collision term, and 
a pure loss reaction equation; the equations are self-consistently coupled and they 
are meant to describe the evolution in time of the phase space densities in a binary 
mixture where the two species interact in such a way that each species inhibits the 
activity of the other, with an interaction proportional to their macroscopic (position 
space) densities. In particular, we want to rigorously derive the system of equations 
in a suitable asymptotics for the particle system. 

The system of equations we shall deal with is a sort of very simple reactive 
transport system, though without conservation of masses; more complex reactive 
transport systems arise in many different contexts, like for instance in the modeling 
of biological systems, of porous media, of radiative transfer or, more in general, of 
various systems in the presence of chemical reactions (see e.g.|SDLj iMMj). We 
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2 CEDRIC BERNARDIN AND VALERIA RICCI 

are interested here in the analysis of the simplest reactive self-consistent coupling 
(i.e. the absorption coupling), and we shall therefore not include other terms in the 
equations. 

From the point of view of particle modeling, an analysis of the absorption self- 
consistent coupling for one species nonlinear equations can be found in models 
for reaction-diffusion equations, such as the ones proposed in |NORj and [Szn] : 
in both papers, the particle system evolves according to a Brownian motion (to 
get the diffusion) and destruction of particles occurs with different mechanism: in 
the first paper the destruction is stochastic, with a death rate for the Brownian 
particles which is a function of stochastic exponential times, in the second one it is 
deterministic, and it happens as soon as collisions between particles occur. 

We shall consider a binary semi-deterministic system where both kinds of inter- 
actions occur. More precisely, the particle system we shall consider consists in point 
like (light) particles moving uniformly among fixed, spherical particles (obstacles). 
Particles not belonging to the same species do not interact among themselves and 
particles of different species interact in the following way: a light particle becomes 
inactive (or is absorbed) at the first time it meets an obstacle and an obstacle 
becomes inactive (or disappears) in a stochastic interval of time whose size is con- 
nected, through a local mean-field type interaction, to the number of light particles 
traveling within the area of detection (range) of the obstacle. 

The main difficulty in the derivation of the (otherwise simple) system of partial 
differential equations originates from the self consistent structure of the problem: 
in order to overcome the mathematical troubles introduced from the self-consistent 
terms, we shall adopt a natural scheme for facing self-consistent problems in partial 
differential equations and particle systems and we shall prove the convergence using 
two levels of approximation of the original particle system, the first one eliminating 
the self-consistent structure and the second one reducing the correlations between 
the two species with respect to the original many-particle system. Our formalism 
is quite explicit and the strategy we shall adopt is similar to the procedure used 
in [NOR] . The adjustment of this procedure to our particle system is not trivial, 
since, at variance with the situation treated in |NORj . we deal with a two species 
system interacting in a non symmetric way and the instantaneous absorption of 
the light particles is in some sense a singular interaction with respect to the local 
mean-field type interaction which govern the obstacles lifetime, i.e. the interaction 
type considered in [NOR] . The presence of the deterministic component (the light 
particles) imposes a sufficiently careful analysis of the trajectories, in the style of 
the analysis performed for particle models of linear equations based on similar de- 
terministic components (see e.g. |BGW| . [DRTJ|DR2]). Nevertheless, the structure 
of the stochastic component (the obstacles) simplifies by a considerable amount the 
mathematics with respect to a two component, totally deterministic system. 

Our final theorem establishes a weak law of large numbers for the empirical 
measures (associated to the two species of particles in the mixture) to the solution of 
the system of partial differential equations, almost surely with respect to the initial 
distribution of positions in the phase space of the light particles and in probability 
with respect to the distributions of positions and life time of the obstacles. This 
weak law of large numbers is valid in an asymptotics where the radius and the 
effective range of the obstacles vanish, keeping a finite action in the limit, and the 
number of particles grows up to infinity, these quantities being related in a way 
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that we shall determine while proving the theorem. As a part of this relationship, 
we shall prove a simple condition (relating the diminishing rate of the effective 
range of the obstacles to the number of light particles) which guarantees the weak 
convergence of the product of an empirical measure times a sufficiently regular 
approximant of a Dirac delta distribution toward the product of their weak limits 
(under suitable regularity assumptions on the weak limit of the empirical measure 
and on the choice of the delta's approximant). This condition entails very useful 
bounds for our estimates and allows to get easily the required asymptotic result. 

2. The equations, the particles model and the main result 

Throughout the paper we shall use the following notations: in dimension d > 2, 
we denote by a; G K'^ the position of a light particle and by w G M'' its velocity; 
t € IR+ is the time variable. In general, configurations of M variables in R^ will be 
denoted by boldface letters with subscript M [jm = (j/i, • • • , Vm)) and sequences 
of variables by capital letters with the subscript oo (Foo ~ {2/i}i^i)- When needed 
(e.g., in the functional spaces) we shall use R^ and R^ resp. for the position and 
the velocity spaces. For p G R'^ and r G R+ we shall denote by Br{p) = {y G R'' : 
\p^y\ < r} the ball of radius r centered in p. For a given z = (x, v) G R'' x R'^, the 
free flow associated to the light particle with initial position (in the phase space) z 
is 

T*(z) = (ri*(z),T2*(z)) = {x + vt,v), t G R+, 

while x{t) = Tl{z), t G R+, is its trajectory, and for e > 

(1) Teit,z)^{yeR'' ;3se[0,t), \y ~ x{s)\ < e} 

denotes the flow tube of radius e associated to the trajectory up to the time t. 

Unless differently stated, for a given random variable r;, we denote by P^ its 
associated probability distribution, and for a measure n (random variable rj), we 
denote by E^ (E,,) the expectation w.r.t. vr {rj). In order to simplify the notation, 
in many stochastic variables depending on configurations ym, we shall label this 
dependence by a simple subscript M, instead of rewriting each time the whole 
configuration. 

We shall moreover denote ^ the weak convergence (convergence in law) in the 
space of finite measures and ^' the *-weak (vague) convergence on the space of 
Radon measures, and by K (with some subscript) any generic constant whose value 
needs not to be specified. We shall sometimes shorten the notation for sums and 
difference of functions with the same argument as (/ ± g){'w) = f(w) ± g{w). 

For r > 0, we consider the following system of partial differential equations for 
the two densities / = f{t, x, w), a = a{t, x), t e [0,T], x e R'' , v e R'' : 



(2) 



dtf + v-V,f = -Cd\v\'^f 
dtu = -e(/jjrf dvf)a 
f{0,x,v) = fo(x,v) 
a{0,x) = ao{x), 



where Q and Cd = //a;PK''-itj|=i> I'^'^^M'^i l*^! — Ij are positive constants. In order to 
have a convenient existence and uniqueness theorem for the solution (cf. Appendix 
13]), we assume /o e L^iRi;W^^°°(Ri)), vfo G L\Ri;W^^°°{Ri))r]L°°iRi x R^J), 
v^fo G L\Ri;L°°{Ri)) and ao G W^-°°(Ri). 
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We shall show that the system ([2]) can be derived from a semi-deterministic 
particle system of the kind we specified in the Introduction. The particle system 
will be described in next paragraphs. 

Obviously, once established the correct asymptotics to get equation ([2]) from the 
particle system we have chosen, we shall simply write lim to denote a limit in 

n— f oo 

this asymptotics (n and e being related). This peculiar use of the notation will be 
consequent to the hypothesis in the environments where the notation is used. 

2.1. The particle system: initial data statistics. We consider configurations 
of spherical fixed obstacles of radius e with stochastically distributed positions at 
time t = and we denote by cm = (ci, . . . , cm) the coordinates of their centers. 
Obstacles may overlap (i.e., configurations such that for some i, fc, \ci — Ck\ < 2e are 
allowed) and M = corresponds to absence of obstacles. In this paper, we shall 
assume that obstacles positions follow a Poisson distribution with parameter /Zg, i.e. 
that the probability distribution of finding M obstacles in a bounded measurable 
set A C M'' is given by: 

M 

(3) PidcM) = e^'^^I^l^ fe^^i • ■ • '^^A'^' 

where (here and in what follows) [A]l denotes the Lebesgue measure of the set 
A. We shall adapt the initial datum cro{x) to this choice for the statistics of the 
obstacles. 

We consider then n point-like particles, located initially at positions xi,. . . ,Xn 
and moving uniformly among the obstacles with velocities vi, . . . ,v„; we shall de- 
note the phase space coordinate of the i-th light particles as Zi = (xi,Vi). A point 
in the n-particles phase space is denoted as z„ = (xi, wi, . . , a:„, «„) and a sequence 
of initial data is denoted as Z^o = {zi}°^^ G [W^ x W^)°° . 

We describe both species of particles by means of their empirical measure (see 
e.g. [Gi]). 

Given a 1-particle probability density /o, we denote by S^ the (infinite product) 
probability measure defined on the space of infinite sequences {W^ x M"^)"" by /o, 
i.e. such that, for n = 1, 2, . . . and A'^^ x A^*^ c M'' x W^, 

n „ 

(4) ^(Z^ : z, e A« X A«, 1 < J < n) = rr / /o(x, v)dxdv. 

The sequence of initial empirical measures for the light particles {m^}^i is then 
such that: 

1 " 

(5) n'^{t,x,v;Zn) ^ -^^5z^{x,v) ^ fo{x,v) ^ - a.e. w.r.t.Zoo- 

z— 1 

We shall require moreover some regularity condition on /o. 

2.2. The particle system: dynamics. We define now the dynamics of the sys- 
tem. 

We consider a sequence Too = {ti, . . . , r/c, . . .} of independently distributed ex- 
ponential variables and we define the following stochastic functions: 
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The risk function at a given position c, 



i-t 



1 r 

(6) Vn,e,M{t,c) = - y] / dsqn{Xi{s) - c)^n,s,M i^' ^i)' 

where qn{x) = af^q{anx), is an (at least continuous) approximant, up to 
a multiplicative constant, of the Dirac delta distribution, with q a non 
negative, radial {q{x) = q{\x\)) function such that J^^ q{v)dy = > and 
in,e,M is the stochastic function defined below. 
• The life functions \j, resp. £,n,e,M{t, z) for a light particle with initial 
position and velocity {x,v) and ??ra,£,i\f (^jCfc) for an obstacle centered in 
Ck e {ci,...,cm}, 

fn,£,o(t, z) = 1 



(7) 



M > 1 

£,n,e.M{t,z) — I M 

{x(s)i U B^{ch)ri^,^,M(s,Ch) Vse[0,t)} 
h=l 

• The maximal collision time between a light particle with initial phase space 
position z = {x, v) and an obstacle located in c, 

(8) T,,, = inf {s: \x{s) - c\ < e}, 

with no reference to the activity of both particles. Because no life functions 
are involved in this definition, Tj c can be infinite. 

In the particle system defined through © and ([7]) , the obstacles become inactive 
at stochastically distributed times defined at a given position Ck through the risk 
function ([6]) as the time tk = tk{ck,Tk) s.t. V{tk,Ck) = Tk, while the light particles 
are absorbed at the first time they meet an active obstacle (i.e. as soon as \xi{Tcou)~ 
Ck\ = £ for some k and TcoU < tk] notice that TcoU is defined w.r.t. active particles, 
and it is therefore a different variable w.r.t. T^^c)- 

Unless we require q to have compact support, a light particle located at x can 
affect the life function of a far obstacle (it suffices that q{x — c) 7^ 0); on the contrary, 
a given light particle interacts only with obstacles which are met by its trajectory 
in space. Since, for a given n, the volume including all light particle trajectories up 
to time T, V„, is such that 

(9) Vn C -B jnax |3;i|+ _max |i;j |T+1 (0) = -Bji, 

and, in studying the interacting system, we do not need to consider quantities 
related to obstacles which can not be met from any light particles, we may (for a 
given n) restrict the expectation value with respect to ^ to the volume A = A„, for 



M 
We use, here and later, the notation |J B^(ci^)r]n^i,^M{s,ci^) to mean (J B^{cf^). 

fc = l ' ' fce{l,...,M}; 

Vn,e,Mi.'>,Ck)^l 

Although not formally correct, this notation allows us to write less cumbersome formulas. 
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a growing sequence of bounded Lebesgue measurable sets {A„}„ s.t. i3„ C A„ and 
lim A„ = K*^. We denote by E" the expectation value with respect to the centers 

n— )-oo 

distribution Q with such a choice for A. We shall denote by E" = EJ?Er = ErE;?, 
where E7- is the expectation value w.r.t. to the exponential times sequence T^o, 
and by P" the corresponding probability distributions. 

2.3. Scaling laws: heuristics. In order to derive the system ([2]) from the particle 
system defined above in the macroscopic limit where the radius e of the obstacles 
vanishes, we have to choose the scaling for densities of both species of particles in 
a suitable way. 

As a first requirement, we want the mean free path of the light particles to be 
finite at the scale at which the system is observed, in such a way to keep track, 
in the chosen asymptotics, of the interaction between the light particles and the 
obstacles. To this purpose, we fix the rate /i^ of the Poisson process ([3]) to be such 
that: 

(10) ti^e'^-^ = M > 

and, since we consider a uniform initial macroscopic distribution of obstacles, in 
(m we assume aQ{x) — fi. 

When Vn,e.M = (i.e. g = 0), formula pO)) defines the so called Boltzmann- 
Grad scaling, which has been analysed for linear particle systems in particular 
in connection with the asymptotics of the Lorentz gas and its variants (see e.g. 
[GallSg lBBS] . [DRl] for the stochastic case, }BGWl[G2llRW] for the periodic case 
and the beautiful review |G3| -focused on the periodic case-). When Vn^e,M ^ 0, 
(fTO| guarantees the finiteness of the mean free path for a given light particle as soon 
as Vn,e,M, in flow tubes and for a non- negligible set of obstacles configurations, 
is bounded uniformly in n. This is true, under the conditions on /o specified 
at the beginning of Section [2] and the conditions on q given in Section 12.21 (i.e. 
q G L'^{W^) n C(R'^)), whenever a^n's < const. 

A second requirement is to have, in the £ — ?■ asymptotics, vanishing correlations 
between light particles and obstacles, so as to obtain a Markovian limit (i.e. without 
memory effects), coherently with the structure of ([2]). Since the dominant part 
of correlations is associated to grazing crossings of light particle trajectories, and 
in particular to the mean volume V^ spanned from grazing trajectories crossings 
within the effective range of an obstacle, this requirement connects the scaling in 
the effective range a~^ to the scaling in the mean density of the obstacles. In order 
to have negligible correlations, the mean number of multiple collisions per unit 
volume has to vanish in the limit. This leads to the condition ji^a'^Vg — )► 0. Notice 
that possibly Vg = 0{e'''~^), and having a better asymptotic behavior depends 
essentially on the regularity of the limit measure of the light particle component. 

We have then a third, technical, requirement due to the form of (|6]): since even- 
tually we want to obtain a mean-field limit, and therefore we want to be able to 
perform limits of mean values with respect to //° of sequences of functions con- 
verging to singular limits (such as Dirac delta distributions), we need a condition 
assuring the convergence (in a suitable sense) of the product of such sequences times 
the initial empirical measure: this is achieved if the empirical measure converges 
faster to its (sufficiently regular) limit than the chosen approximant concentrates, 
in such a way that the empirical measure is in practice equivalent to its limit den- 
sity well before than the delta's approximant concentrates in its center. We can 
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guess roughly that this happens if the fraction of particles fluctuating around mean 
values in a volume corresponding to the effective range of an obstacle vanishes in 
the chosen asymptotics, and therefore -j= — o{a^'^). We shall prove a more precise 
asymptotics in Appendix 16.21 

Given ([T0| and the just described scaling laws, for a given configuration of obsta- 
cles, the empirical measure at time t for the light particle component (representing 
its mesoscopic density in phase space) is: 

1 " 

(11) Mn(^, a;, W;Z„,Cm) = -y^ STt(zi]{x,v)^n.s,M{t, Zi) 

and the (macroscopic) density of the obstacle component (i.e. the number of ob- 
stacles per unit volume) is expressed as: 

M 

(12) Cr„(i,X;Z„,CM) = £''~^y^^efc(a:)?7n,e,Af(t,Cfc). 

k=l 

2.4. The limit process. Under the conditions we assumed on (/cctq), as shown 
in paragraph l6.1l in the Appendix, the problem ^^ admits a unique solution (/, g) £ 
L°°([0, T]xR'^x M'^) X i°°([0, T] x R'*), with / dvf G i°°([0, T] x M'') (and actually 
J dvf{t, -,«) G Cb(M'*)). This solution can be expressed in semi-explicit form as: 

fit, X, V) ^ foix- Vt, W)e-'^''l^l -fo dsais,x~v{t-s)) 

(13) 

ait,x) = ao(2;)e-®^" '^"^"''''"/("■^'"^ 

In order to be able to compare solutions of ^ with the stochastic measures ITTI) 
and (J12p . we define the risk functions and the life functions associated to the limit 
process described by ^ as: 

V/(i,c) = e f ds f dvf{s,c,v) 
(14) ^° ^^' 

Vi{t,c) ^ I{Vi^(t,c)<T} 

Ul{t,z)^C,\Ti{z)\ f dsa{s,Tl{z)) 



(15) 

where r and Tp are independently distributed exponential variables and r = t^ 
when c = Ck- 

In this way, the semi-explicit form of the solution (J13p can be re-expressed as 

f{t,x,v) = foix - vt,v)ErM{t,T-\x,v))] 

(16) 

a{t,x) = (To{x)Er[iil{t,x)]. 

In order to simplify the notation, we shall omit in what follows the dependence 
on the density functions (Vl — VI,Ul — U1). 
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We want to establish a (weak) law of large numbers for the measures pT|) and 
([T2I1 . More precisely, we shall prove the following theorem (we recall that J^iW^) 
denotes the space of C°°(R") functions of rapid decay at infinity [SchllDL] ): 

Theorem 1. Consider the non-negative functions /q and q and assume 

• /o G yO^'^ X m.'^)C\L^{Ri;W^-°°{Ri)) IS a probability density such that 

(17) / di;/o e ^(M'') with I dvfo{0,v)>0, 



andvfo e Li(R^;W^i'°°(R^)), v\fo G ii(M^;i°°(R^)); 

q is a radial function s.t. q G S^{M.'^) and J„^ dxq{x) — Q > Q; 

{flnl^i is such that an > 0, lim a„ = cx) and there exists some k g (0, ^) 



such that 

(18) lim 4 = 0; 

n— >-oo 77, 

• {£}^=i = {£n}n=i s-i- £n > and 

(19) lim a^e'^ = 

n— >-oo 

for some C £ (0, i - ^). 

Then, ^-almost everywhere w.r.t. sequences of initial data Z ao and in probability 
w.r.t. P and Pr, when n — ?• oo 

1 " 

(20) fln{t,X,V;Zn, Cm) = -y^STt(zA{x,v)^n,e.M{t,Zi) -^ f{t,X,v) 

n ^^ 

i—l 

M 

(21) a„{t, x;Zn, cm) = e'*"^ ^ (5c,(a;)77„,e,M(t,Ci) ^ cr{t,x), 

where ^„ e M; Vn e m cire defined in varaaravh \2.2\ and {f,cr) is the unique solution 
of 

' dtf + v-\/,f = -Cd\v\'^f 

(nn\ J dt(y = -©(/«<< dvf)(7 

^ ' I f{0,x,v) = fo{x,v) 

^ a{0,x) ^ fi. 

We shall see in paragraph l6.2l in the appendix that hypothesis (fT8)) . together with 
suitable regularity assumptions on /o (all included in our theorem) , guarantees the 
following ^-a.s. convergence 

(23) |«|V°(a;,w)g„(x) - Q\v\^foix,v)So{x) 

blV° -^ \v\'fo 

for j = 0, 1, 2, . . . , F, fc = 1, . . . Q, with given positive integers P, Q. 
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3. Definitions of the approximating systems and proof of the main 

THEOREM 

The main difficulty in studying ^ and its associated particle system ([B]),© is 
that we have to deal with a self-consistent problem. The ffist step is therefore to 
find, both for the limit system ([2]) and for the particle system defined from (|6]),([7]), 
suitable approximating systems which do not share this self-consistent structure. 

We first recall here without proof, since it will be useful for the sequel and we 
shall use it largely in various steps of the convergence proof, lemma 3.2 in j NOR] , 
concerning bounds of distances of stochastic variables of the form ri{t) — I{s(t)<T} 
in terms of distances between their associated risk function S. 

Lemma 1. Let t > be an exponentially distributed real variable and define, for 
i — 1,2, t £ [0,T], rii{t) = ItSitt)<T}: where Si,S are non negative, non decreasing, 
right continuous random functions. Then, for any 5 > Q 

(24) hiW - r;2(i)l < \\Sl{t) ~ S{t)\ + i|52(t) - S{t)\ + /{|5(t)-r|<5}. 

Remark 1. Whenever at least one or both risk functions Si, i — 1,2 are in- 
dependent of the exponential time t, the bound simplifies in the first case as 
^r\vi{t) ~ V^it)] < jEr\Si{t) — S2{t)\ + 26, and in the second case as Er\rii{t) — 
V2{t)\ < \Siit) ~ S2{t)\ . 

3.1. Approximation of the limit system. As shown in Appendix 16. 1[ under 
suitable hypothesis on /q, the solution of the system of equations ([2]) can be obtained 
as the fc — > oo limit of the sequence of solutions of the sequence of linear systems 
defined as: 

(25) f<^'^\t,x,v)^foix-vt,v), a(")(i,x)=/x 

^.^^^-(e^.d.;/^^-!))^^ ^^^ 

fiO,x,v) = fo{x,v) 
^ a{0,x) = fi 

More precisely, the sequence of semiexplicit solutions of ((25|) is, for fc = 1, 2, . . ., 

/W(t, X, V) = /o(x - Vt, v)e-''^\"\ /o dsai^-^Hs^^-vit-s)) 

(26) 

af'^) (i, x) = Me-^ ^0 ds /«. d./c--!) (s,.,v) ^ 

and we have, when fc — )■ oo (cf. Appendix 16. II) 

11/ - f^''' \\L°°([0,T]xR'ixR'i) -^0, II /jjd dv{f - /( '')\\l°°([0,T]xR'^) ^ O7 

(27) 



kJ — (T 



\\l°°{[o,t]xr<') -^ 0- 



The risk functions associated to the fc-th system p5|) are defined as 



(28) V^''\t,c) = e ds dvf^^-^\s,c,v) 



(29) t/«(i,z) = Cd\n{z)\ / dsa(^-^\s,Tt{z)), 

Jo 
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and, because we shall need it later (cf. the definition of the system ([33l) ). we define 
F*^°-' (t, c) = 0. Their associated life functions are: 

(30) 

with T and Tp exponentially distributed times. Of course, (1261) can be expressed in 
terms of V^^'^ and t/*^'^' in a form analogous to (|16p . 



3.2. Approximation of the particle system. In the same spirit, we may also 
approximate, for each n, the particle system described by (jH]), © by a suitable 
sequence of systems. For an initial datum z„ for the n particles phase space position 
and configuration of obstacles cm {M — in the absence of obstacles) , we define 
this system, for integers M, k > 1, j ^ 1, . . . ,n and £ = 1, . . . , M, as: 

K!J,M(i> ci) = 0, vl%{t, q) = 1, il%{t, z,) = 1 

1 ji 

(31) K!.Sm(*' C£) = -g /o dsqn{x^{s) - Ci)^i^;)^i{s, Zi) 

fC") u z) = I 

This sequence is a linearization of the original particle system (which is its formal 
limit when k — >■ oo) in the same way as psp is a linearization of ([2]) and it satisfies 
what is called in |NORj sandwiching property (see ()48p ). This property implies, in 
a given asymptotics for n and e: 

E" 1 Ki^i M - ■^^''^ I ^ => E" I K,e,M - Fl I ^ 



(32) 



E" lC,M - ^"^'^ I ^ =^ E"E,^ |e„,e.M - a I ^ 

E" I / dz^,^^{S!^^J^J - ^C^) ) K =^ E"E,^ \Jdzfin<l>{Cn.e,M - ^l) H 

and using these implications we shall be able to bypass the direct evaluation of 
quantities related to the particle system ([5]), ([7]). 

The advantage in dealing with the two approximating sequences of systems de- 
fined by (PS)) and (I5T|) instead of the original systems is that, for each k, the two 
components evolve in a given field, associated to the functions defined at the pre- 
vious step fc — 1 in the sequence, and the original self-consistent structure is lost. 
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The next step would be then to show that (|25p and (|3T|) are, for a given k, 
asymptoticaUy equivalent. Unfortunately, the system defined from pi|) . when e 
is positive, still keeps strong correlations between all light particles and obstacles 
positions (in the phase space), correlations which are absent in (|25p: this makes 
hard the direct comparison of the two systems. Therefore, we need to define an 
intermediate system in which correlations among light particles and obstacles are 
further reduced. We shall prove then that this system is equivalent in the limit 
n -> oo both to ^5^ and to (pij) . 

3.3. A system asymptotically equivalent to both approximating systems. 

For an initial datum z„ for the n particles phase space position, and for configura- 
tions of obstacles cm [M = in the absence of obstacles), we define, for integers 
^I, k > I, j — 1, ■ ■ ■ ,n and £ = 1, . . . ,M, the intermediate system in the following 
way 



(33) 



in,e,M\'^y'^j) ^ \ 



{^j(s)^^U^B.(ch)/^^(fc„i)(^_^^,^^^j VsG[0,t)} 



ni=i 



^nl,M^^' C^) = - j^ /p ds qn{x^{s) - Q)li^^ jl(s, Z,) 



.(fc) 



where Tk, 1 < k < M, are independent exponentially distributed times. 

In ([33|). the life function of a light particle at level k is defined through fictitious 
obstacles life functions (^/\/(fc-i)(t c)<Tc})i corresponding to the obstacles life func- 
tions at level fc — 1 associated to (l^5|) : the life function of an obstacle at level k, 
^i aW ^ \: is defined through the light particles life functions at level k — 1. In 

this way, the correlation between light particles and obstacles is weaker, compared 
to the same correlation in system pip , and this allows us to prove, when n and e 
verify conditions ^TE\\ and (|19p . the convergence to both systems (P5|) and ([51]) in 
quadratic mean w.r.t. the expectation value E". 

3.4. Proof of the main theorem. We proceed now with the proof of our main 
theorem. To this purpose, we assume we have already proved the sandwiching 
property for the system ([51]) and the asymptotic equivalence of this last system 
to (^5]) : we shall postpone to next sections the proof of lemmas and propositions 
concerning these two facts, since they are the core of the derivation; as pointed out 
in section [3TH the convergence of ([5T|) to (P5|) will be obtained passing through the 
equivalence in this asymptotics of ((33|) to both systems (1311) and (pS]) . 



Proof of Theorem]^ Given the definitions in the previous paragraphs, we write 
and 
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Then, using formula (HS]) for / and ^ for /C^) , we can write, for all (/) G Cb(E'* xE"*) 
and for all fc > 1: 



E" 



(34) 



(35) 



1 " f 

- V'(?!)(T*(zi))^„.e_M(i, Zj) - / dxdv(f){x,v)f{t,x,v)\ 

n 

-^0(T*(z,))(w,M-a)(t,2O| 



< 



E"E^ 



\4>\\oo { E"E,^ 



^Ei^--^"^'^^ !(*'-.) 



i=i 



dxdvCd\v\fQ{x,v) / ds|cr-o-('' ^^ |(s, a;(s)) I + 



(36) 



E"E^ 



1 " _ r 

I - V (/.(T*(z,))e*'^ (t, z,) -/ dxdt^/C^-) (i, X, t;)0(x, «) I 



The term (p4|) vanishes on a set of full measure w.r.t. ^ as a consequence 
of the sandwiching property valid for ((3T|) (cf. Corollary [l] Sec. 14. 3p and of the 
convergence of (I5T|) to ([25]) in the asymptotics specified in (IT51) , (fT91) (cf . Proposition 
a Sec. El). 

Thanks to (P51) . and in particular to the convergence l^l/i^ ^ I'^^l/oj the term 
in curly brackets ([55)1 is bounded (excepted on a set of zero measure with respect 
to 3^) by -K^allcr — (y'^^~^^\\L°°([o,T\y^9.<^)^ SO that, by choosing a suitable fc, it can be 
made arbitrarily small because of ([ST)) . 

Using Cauchy-Schwarz's inequality, we get then for (|5S]) 

2 



E"E^ 



1 _ /■ 

I - V (/.(T*(z,))e('=) (t, z.) -/ dxdvf'^''^ (i, X, «),^(x, v) I 



< 



1 " 



+ / dxdvf^''\t,x,v)(f>{x,v)x 



dxdvf^'^\t,X,v)ct>{x,v) - - V0(r*(z,))IErJ^"('=n^,^O] 



Because of the hypothesis on /q, (/)(T*(-))E^J^('=)(i, •)] e C6(M'' x M''), so that, 
thanks to ([^ (and in particular to /i^ ^ /o), on a set of full measure w.r.t. ^ we 
have 



1 

n ^ — ^ 



-^ 



j=l 



dxdvf^ ' (t, a;, u)(/)(x, u). 



and ((55)) vanishes on this set. 
In the same way, by writing 



Vn,s,M = iVn,e,M ^ Vl) + {VL - v'^'"^) + V^'"'' 
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and 



rW^ 



a = {a — a^ ' ) + a 
we get, for all ip G Ck(]R'') and for all A; > 1, 

M 



(k) 






< 



(37) 



(38) 



(39) 



E" 



M 



e'^ ^y^Jsupp4,{ct)\rin,s,M -riL\it,Ct) 



M 



e'^-^E^-PP'A(^»)l'?^'^-^^l(^'^») 



+ / da;|cr('=) -(T\{t,x) 

•^SUpp-i/j y 



+E" 



M „ 

£'^-iVV'(c,)^^'=)(i,c,)- / dxij{x)a^''Ht,x)\ 



where the term in round brackets psp is (everywhere) bounded by ii'(,|| J dv{f — 
/''^"^•')llL°°([o,T]xR<i) and can be made arbitrarily small because of ([27]). 

The term (157)) vanishes on a set of full measure w.r.t. ^ for the same reasons 
as ([M]) . i.e. as a consequence of the sandwiching property valid for ([51]) and of the 
convergence of pil) to (l?5t in the asymptotics specified in (fTS)) . (fTO|) . 

As for ([M]) . we can write 

2 



E" 



A/ 

e'^-i V V(c.)?7^'^ (i, q) - / dxi^ix)a(''^ (t, x) I 



< 



dc^{c)<7^^\t,c)~ / dcV'(c)CT('=)(t,c) 



-d-l 



dcV''(c)a(*^(i,c). 



Since V e Ca:(M''), we have ?A''(-)o-('=)(i, ■) e ^^(K'^)) for /i = 1, 2 and when n grows 
to infinity we have both A„ ^ W^ and e — > 0, so that this term vanishes, in the 
chosen asymptotics. 

Collecting all the assertions about the different terms, the theorem is proved. D 

4. Asymptotic equivalence of ([5T|) and (P5|) and convergence of the 

PARTICLE SYSTEM TO THE LIMIT SYSTEM 

We collect in the present section the Lemmas and Propositions which will help 
us to build up the last step of the proof of our main theorem. 

4.1. A few useful bounds. Let us list a few bounds and formulas which we shall 
use often in our calculations. 

It will be useful to adopt in next sections, for stochastic variables \^^-^' of the 
form 

A "' (t, Zj) ^ I M , 
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(with J-n,M{s,Ch) = -^CAj,TAf,z„(s,c/i) a given random function) the foUowing rep- 
resentation: 

M 

(40) A^--(i,z,) - n (l - Ac.erat...)}^{^„,M(T.,,.,,c.)<r.,}, 

h=l 

-p{l) -p{2) i-rm, 

Given any couple of stochastic variables A "•"^, A ""^ of the form ([301), using 
the trivial inequality 

\X{1a,^\{IbA<Y.\Ia, - Ib,1 

h h h 

for M > 1, we get the bound 

M 



(41) 
where 



^1 ^ -^2 



{t,zj) <J2^^'^^r,it.z,)\v^--'' - V^"'''\iT.^,c,,Ch), 



h=l 



r?^--(r.,,c„c,) = Vil(T.,,.,,c.)<..j' * = 1'2. 



The bound (|4T|) will be largely used in the last paragraph of this section. 

Moreover, because of the hypothesis on /q (assigning initially uniformly bounded 
number and kinetic energy limit densities for the light particle component, cf. Re- 
mark [HI in Appendix 16.21) . ^-a.e. with respect to Z^o, we may use the bounds 



1 " 

(42) l|-E9"(^J-(^)-^)ll- ^ ^1 

"^3 = 1 
1 " 

(43) \\-Y,qn{xj{s)-c)[re{s,Zj)]L\\oo < e''-^CdTK2 

where Ki,K2 > are constant independent of n. 

4.2. Quantities related to flow tubes. When bounding correlations in our parti- 
cle system, we shall need to evaluate expectation values with respect to the Poisson 
distribution ^ and to /i° on volumes corresponding to the intersections of flow 
tubes. 

Defining, for two vectors v,w G M'', cosa(i',M;) = -r^Ti^ (and a^ = a{vi,Vj)), 
and recalling the notation Br{p) introduced at the beginning of Section [2] to denote 
the ball of radius r centered in p, we have, for the intersection of two flow tubes 
associated to the particles with initial phase space positions Zi and Zj , the following 
trivial bound: 

Lemma 2. For any /3 e (0, ^).' 

(44) [re{t,Z,)nTe{t,Z,)]L< 

[S2(0)]Le''^'-^)/|sina,,|>e^+C^de'~'min(|i;,|,|i;,|)r/|,i,„^^l<,,. 
Proof. We have always 

[Te{t,Zj)nre{t,Z,)]L < inm{[Te{t,Z,)]L,[Te{t,Z,)]L) 
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and, whenever sinaij ^ 0, denoting as y the (unique) crossing point for the trajec- 
tories Xi{s), Xj{s) (i.e. y = Xi + ViSi — Xj + VjS2 for some si, S2 € M), 

Teit,Zj)nTe{t,z,) cB^^iy). 
We can therefore write: 



(45) \%it,Zj)n%it,z,)\<\B2m(-^) ACds''~'mmi\v,l\v,\)t 

\ sm dij J 

and (gH) fohows. 

D 



Thanks to the weak convergence of the initial empirical measure toward a regular 
density, we can estimate the measure (w.r.t. /^^ ) of the set corresponding, for a 
given velocity w G M.'^, to the grazing collisions. We can prove in fact the following 
Lemma: 

Lemma 3. Let /i^ ^ /o € y{M.'^ x M'*), l3 > and w € R'^. Assume H^) . 
Denoting cosa^ = J^TUTu]' ^^^ following bound is verified 



1 " 

-^I\sina,\<ef< < Ks^ + o{a:^/) 

foranyCe{0,^P). 



d-l , 



n 

i=l 



Proof. We observe first that we have lim "''"^'"'^ I^^^jIj = 0, since /o G. 

n— >oo " 

so that, uniformly in e, the contribution of the corresponding term vanishes when 
n -^ oo and it is actually o{a~'^) in the prescribed asymptotic^- 

We now evaluate the contribution of grazing crossings of particle trajectories. 
We choose a suitable (standard) C°° regularization R^ of /|„.^|^|„||^|, for instance 
R^ is s.t.: 



»(",»") 



Here the parameter 5 denotes the radius of the set where the regularization differs 
from the original characteristic functions (i.e. the set in W^ x W^ s.t. | sina(u, w)\ < 
(5), and, for all (5 > and < t < ^^, we write: 



1 v-^ 1 V \ £*■" 



-./3 / ^.^., U^^^) , r^_ ^"^ 



< e'P / dxdv -"^ ' ' + Kr-— + o{5). 

^ \v\'^\w\'^ ' 

The error coming from the regularization of the characteristic function /u.iu|^u||iu| 
is o[S) because of the weak convergence of /i° toward /o G ^(M'^ x R'') and, for the 
same reason, all constants (here and in the following estimates) are uniform in w 



\n fact, since by Toeplitz's lemma, for all a G (0,1], - Y^ ^ #{J>j-"'l ['"ll^jl) -^a vanishes 



1 1.1, #ljH^7 "^1 — 1^1 1^7 I } / 1 \ \ 

when n — > cxd, we have at least — = o{^=) ). 



k-1 
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and n and depend only on the dimension d and on few L^ norms of /q. We obtain 
then, for suitable choices of S, 



1 " 

-2_^^\yi-w\T^\vi\\w\I\sinai\<el^ - -^^ ' 



n 

i=l 



for any ( G (0, t/3) and therefore for any C e (0, %^/3)- Collecting the two estimates, 
the lemma is proved. D 

We shah use Lemma |3] for free choices of /3 e (0, ^), so that ( G (0, ^P) C 

Remark 2. The estimate obtained in Lemma |3] is valid for any given sequence of 
empirical measures {^°} such that /x" ^- /o G y{R'^ x M"^) and it is actually a 
stronger assertion with respect to what we need to prove Theorem [1] (which is 
valid ^ — a.e. w.r.t. Zoo)- In order to prove Lemma |3l where, as stated in the 
introductory sentence to the lemma, we bound the measure of a set with respect 
to the empirical measure /i°, we need to evaluate the measure of the set {{x,v) : 
\v ■ w\ — \w\\v\} with respect to the measure /i° because the measure of this set 
vanishes only asymptotically (i.e. with respect to the limit measure with density 
/o). We choose to evaluate this measure with respect to the parameter afj (i.e. as 
o{a^'^y) for further convenience. 

A similar, simpler statement can be proved ,^-a.e. w.r.t. Z^o'- in this case we 
obtain the bound ^Y^^=i^\sinai\<ei^ ^ Ke'^ (^-a.e. w.r.t. Zoo). This alternative 
statement could be used instead of Lemma [3] to prove our main theorem. We prefer 
nevertheless to use Lemma |3l getting in this way bounds which are (as much as 
possible) valid on the whole subset of initial sequences Zoo such that the associated 
sequence of empirical measures {/x° In^i converges weakly to /o G S^{W^ x W^). 

Remark 3. Notice that, since -^ ^ 1 = o{a^'^) and all bounds are 

i-J-\t>i-Vj\ = \vi\\vj\ 

uniform in the velocity w, we get also 

-. n ^ n ^ n 

ij = l j=l i=l 



A last very useful tool will be the following parametrization of the points c e Te{t,z). 
Let u = Tz.c be the maximal collision time defined in ((S|; then we may define the change 
of variables 

c = x{u) + suj, ue[0,t), uj&dBi{0). 

Using this parametrization, for each non negative function V G I/°°([0, T]; Vy^'°°(R'*)) we 
can write V(Tz,c, c) = V (u, x{u) + elS), getting then 

(47) / dce-^^^^-'^''^' = e^-^CM ( f d^[e-^("'"'"» + 0,(m, x{u))]) 

where we may bound the remainder through the trivial inequality |e~^ — e~^|<|a; — ?y|, 
for x,2/ > 0, and the function ^^ is such that supro^j HiJ^ellioojud) = 0(e) — >■ when e — >• 0. 
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4.3. The sandwiching property. We prove now the sandwiching property for the sys- 
tem defined by pi|l and its imphcation on the convergence of ((Tjl. 

Lemma 4. (sandwiching property) Consider Ki,e,M, £,n,e,Ai, Vn.e.M defined by (6j), f^ 
and Ci';l,M, K'Sm, vI^Im defined by {21- Then, fork = 1,2..., 

e{2fc-l) ^ A2k+l) , ^ , A2k) ^ ^(2fe-2) 

'in,s,M — 'in,e,M — ?",e,M r^ Sn.e.Af li ?n,£,Af 

. .o\ T/(2fc-2) , -,.(2fe) ^ , . ^ ^A■2k+l) ^ ■,^(2fe-l) 

(2fe-l) , {2fc+l) ^ , {2fc) , {2fc-2) 

Vn,e,M < %,e,M < Vn,e,M < V„,,,M < Vn,,,M ■ 

Proof. From definitions ©, dZ} and ^, we have = VJJ'i^M < 'Ki,£,m < K!,^i,M. which 
implies f?!"^ j,^ > Vn.,i^,M > Jyi^eM' which in turn implies ^i^^ m < ?",£,*/ < Ci^' ^i < 

We have therefore, using again these definitions: 

?n,£,Af — 'in,£,M — ?»i,e,AI r^ Sn.e.Af ^i Sn,e,M 
(1) ^ (3) ^ ^ (2) ^ (0) 

and (|48)) is valid for fc = 1. 

Since the following chain is also valid: 

(2*:-2) ^ {2fc) ^ ^ (2fc+l) ^ (2fe-l) 

<,e,Af > %,.,Af > '7n,= ,Af > '7„,,,a/ > '7„,,,A/ ^ 

^(2fc-l) t(2fc+l) <^ t <r f(2fe+2) , t(2*:) , 

?n,e,A/ — ?n,e,Af — ?n,e,Af r^ ?n,£,Af — ?n,e,Af ^ 

T/(2fe) < T/(2'=+2) < T/ < T/^^fc+S) < T/(2'=+l) 

"^n.e.A/ ^ ''n,e,A/ ^ Vn,s,M S »'„,£, Af ^ ^n,6,M 

the statement of the proposition follows by induction. D 

We may then prove the following corollary to Lemma |4l 

Corollary 1. Consider V„,e,M, ^n.e,M, r)n,e,M and V„^J}^,i, Ci'^'^Af Vn,l,M defined as in 
Lemma[^ and Vl and ^l given by p^ and 1115}) . with (/, a) unique solution of (^. Then, 
when n — >■ oo,e ^ 0, for any sub-linear operators L on L°°([0,r] x R"^ x R"^) 

MV1':Im - ^*''l ^ ^ L\V„,e,M -Vl\^0 
L'^'r\rin.l.M - V^''^] ^ =^ LEr\Vn,e,M - Vl\ ^ 
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Proof. We denote here by (/>+ = (^ A 0, <^- = —(f) V resp. the positive and the negative 
part of cj>. 

Let Tn,E,M be one among the non negative functions Ki,e,M, 1 — i,n,s,M^ 1 — ??n,e,M, 
J dz^i'^(l)+ ^„,e,M, / dz^j,'^(l>- ^n,£,M and J^n.l.M its approximant of order k (defined through 
l|3ip V Let Tl be its associate limit function (defined through (|f4|) and (jfSp ) and J^'*' the 
approximant of J"l (defined through (f28|) . (f29l) and pOl) ). 

Because of (|48[l . we have, for a given fc > 1: 

T-(2fc) Tr(2'= + 1) TT 

and therefore, 

I Tr(2fc+l) _ Tr(2fc+i) j , I Tr(2fc) _ -c(2k) , 

Since ||Vt-\/'''^||ioo([o,T]x»'i) and ||Erp[^L-^<''^]||ioo([o,T]xK'ixR<') vanish in the fc -^ oo 
limit (cf. ((27l)), and since ^Tp[\riL - r?'''''!] < |Vl - l^'*'] (cf. Remark[l]), the assertion is 
proved. D 

As already pointed out. Corollary [T] is valid in particular when L is an expectation 
value operator. 



4.4. Equivalence between (f33]) and (p5]). 

4.4.1. Motion of light particles in a decaying medium: equivalence of the light particle 
components of {23) and \25\) . We prove here a lemma which describes the behavior of the 
light particle component of the system (|33p in the n — > oo limit. In this lemma, we consider 
a particle system such that the rate of death of the obstacles is given and independent from 
the light particles component, while the light particles are instantaneously absorbed each 
time they collide with an obstacle. The correspondent limit system consists of equations 
which are coupled only in a very weak way and the parameters determining the asymptotics 
are independent. 

Letp(f,c) G L^([0,r]; W^'°°(R'')) be a non negative function. The risk function ^(i, c) 
defined, for t < T, by 

V{t,c) = / g{s,c)ds 



is non decreasing as a function of t. 

As in the previous sections, we define: 

V (t,Ch) = -^{V(t,Cfc)<T,J 



(49) 



Cn,e,M(ii Zi) — I M 

{x,{s)(t: U Se(cfc)^^(s,Cfc) ve[o,t)} 
fc = l 



We have then: 



Lemma 5. Letg{t,c) e L^([0, T]; H/^'°°(R'')) and^Z,e,M be defined by g^. Consider the 
sequence Zoo G (R'' x R'')°° and its associated sequence of empirical measures {/^Jlj^JLi. 
Assume /i° -^ /o G .y and \v\^°„ -^ \v\fo G =5^. Then V(/> G CtiM.'' x 



d s.y TU)d\ 



lim E" 

n—^oo 
e-»0 



1 " f 

-y^C,E,Mi''^^ ^■^)'t>{T\zi)) -/ dxdvf{t,x,v)<j){x,v) 
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(independently of the order of the limits), where (/, a) is the (unique) solution of the 
following system of partial differential equations: 

{ dtf + vV:,f = -\v\Cdaf 
dtu = -g(T 
f{0,x,v) = fo(x,v) 
a{0,x) — fj,. 



Remark 4. Of course, when V — V^'' ^ , we have ^^^ m ~ i„ i M' ^'^ that obviously (by 
the triangular inequality) Leninia[S]implies, when e — >■ 0, n — )■ cx), the limit E"[| J dxdvfj!^{^. 

Proof. According to the definition of maximal collision time ([8]), we have 

V/i : 1 < ft < M, 



u(k) 

n.£.M 



Cn,e,Af (i, 2j) — 1 



T,^,c, At> inf.gM^ {r;^(s,c,,) = 0} A i. 



The function F{-, c) — inf{z : z — V ^(•, c)} is well defined and 

inf{s e R+; r?^(s, ch) = 0} = F{T^^,Ch). 

Moreover, the maximal collision time Tzj,c is greater than t if and only if c does not belong 
to the tube 71 (i, Zj). It follows that 

P"[C,M(i'^j) = ^] = P"[{l<h<M:n^,,,At>Fir,„c,,)At}] 



= E" 



IVI 



e%{t,zjy, T,^.c;.>v-i 



(Te^,Cfc)At}j 



exp 



dee J' 



Te{t,Zj) 



Using (I47p . because of the weak convergence of /^^ and \v\ii'^, we get 



(50) E" 



1 " 

- y~'?r,e,M(i, 20</'(r*(2:»)) 



n— >CXD.£— >0 



B'^xM'' 



dxdvf{t,x, v)<j){x, v) 



the result being independent of the ordering of the limits. 
From (jSUp . we prove Lemma [S] as follows. We call 



fm = E" 



1 " 

— y" cr,e,M(i, 2i)?X£,M(i,2j)<?!'(r*(2:»))0(r'(2j)) 






E^n ~ / dxdvf(t,x,v)(j){x,v)—]S," 
We write then 

E 



- y~'?r,e,M(i, 2i)</'(r*(2i)) 

77 Z ^ 



E" 



-'^^n,EM{t,z^)(j){T^{zi))~ / da;di;/(t, 2;, -y)(/>(a;, w) 

— '^(,n,s,Mit,Zi)(j>'^{T^{Zi)) +fl„+£2n/ dxdv f {t,X,v)(t){x ,v) 
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;^EeL.M(i,^>)'^'(r*(20) 



where E" 

to CT . 

£2n I dxdvf{t,x,v)<f){x,v) 

As for £in we observe that, for i 7^ j: 



when n — )■ 00, obviously vanishes and, thanks 



dxdvf{t, X, v)4>(x, V 



kR'' 



exp 



-/i£ / dee 



-V(T,.a,c) 



Me 



/ 



dee ^ ^3.=' ' 



exp 






-V(Tj 



.,..) 



i.j — 1 



We write 1 = /|sinai,i>e'3 + ^| sma^^Ke^. for /3 < H- Then, using ^ and [7;(t,%) n 
Te{t,Zi))\L < Kfj,£'^~^ (always valid for intersections of flow tubes, as shown also in the 
proof of Lemma [2]) we have: 

1 ^ 

■ ;=i 

'^ — I\s\-nai\<ef<)^ from (fl5)) and Lemma 
1 " 



As for the remaining part (since ^uina l>e 
|2]we get: 



ij=i 



1 " 
o(a-'') +0(e«) + - ^ ,^(T*(zO)0(T'(^,))x 



!.J = 1 



Msf/ dee-^(^^.-=>+/ dee- 



xexp -/is ( / dee"' ^'^' •'=■'-'+/ dee"'^<^^J''='"M+^te^"'"' 

since in this case /iE[7^(i, Zj) P 7^(f, Zi))]i = 0{e^^'^^) — > 0, so that (using (|47)) once more) 

2 



fit 



M'^xR'' 



dxdvf(s, X, v)(j){x, v) 



and Lemma [5] is proved. 



D 



We emphasize again that all results in this Lemma are independent of the ordering of 
the limits n — >■ 00, e — >■ 0. 



PARTICLE MODEL FOR COUPLED PDE'S 21 

4.4.2. Equivalence of the obstacle components of Ii33\} and Ii25\) . Lemma [5] describes the 
asymptotic average behavior of the light particle component of the system (I33|l . The 
asymptotic average behavior of its risk function, for an obstacle located in c, should 
be equivalent to the behavior associated to the risk (|28|l . This fact is described by the 
proposition proved in this section and its corollary (notice that the risk function of a given 
obstacle with position Ch, differently from the generic risk function associated to a given 
position, is defined only if there exists an obstacle in c^ ). The second proposition will be 
useful in the proof of equivalence of (|33|l to (|31|l . in the next section. 

It will be convenient to define the following stochastic functions. For z = (x, v) € 

(51) pf^{i, C, 2) = 1 - /{cere(t,B)} ■'"{V('»)(T,,,,c)<re} 

p {t,z)— du \v ■ uj\dujexp i—V {u,x{u))] 

Jo JaBi(o) ^ ' 

Qn{s,t,c) = g„(xi(s) - c)qn{xj{t) -c) 
K' "-■) (s, i, c) = E., [pf ) {s,c, ^Op^'' {t, c, z, )] . 

Given these definitions, a fundamental tool in all estimates in this section is the fol- 
lowing lemma: 

Lemma 6. Take a sequence Zoo G (R'* x R'')"" and its associated sequence of empirical 
measures {/ij^j^i. Take a sequence {A„}5^i of hounded, Lebesgue measurable sets such 
that Bn C A„ (with Bn defined m l^) and assume /x^ ^ /o G ^ and \v\^q -^ l^^l/o G ^• 
For < s < t, we have 

J. dcTZ7^'\t,s,c) ^l-I^[p(^){t,z,)+p(''\s,z,)il-S,,)] 



^ ' R?\z.,z,,t,s) 

^ [An]L 

where, for each C G (0, ^ — j^), 

1 " 

(53) -^ ||-Rf^(2i,2j)IU-([o,T]x[o,T]) < Kri£''~'^'^'^\vj\ + e'^~'^[o{a~'')] 

i = l 

and 

1 " 

(54) lime^"'^--2 ^ ||7?^^>(2i, 2:j)||ioo([o,T]x[o,T]) = 0. 

Proof. Formula H52p follows, through simple calculations, from H47|) . expanding the expec- 
tation value in the definition of T^.^^ , given in (|51l) . and writing 



-V'<'">(T,,c,c) _ d-1 



p^'''(t,z) + Cd\v\ du(j)e{u,x{u)) 

Jo 



dee 

doing this, we obtain for the remainder 

\R?Hzi,Zj,t,s)\ < e''~^Cd{\vi\ + 1i'j1)|1Ve||loo([o,t1xR'') 
(55) 

+ [%{s,ZJ)nTeit,Z^)]L 

where |lv'ellL°=([o,T]xH'i) = C'(e) -^ 0. 
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We use then the bound (I55p to obtain (|53p and (|54p . In the single sum in (|53p . so 
as in the double sum in (|54p . thanks to the weak convergence of |u|/iJl, the contribution 
of the first term on the right-hand side of (|55p is 0{e''). As for the contribution of 
[Te{s, Zj)r\T'e{t, 2i)]i, we bound the Lebesgue measure of the intersection of the tubes using 
Lemma (2] then, thanks to the weak convergence of /i^ toward the (sufficiently) regular 
function /o, we can use Lemma [3] and (|46|l . with /3 < ^, to estimate the contribution 
coming from the grazing collisions (being the remaining part in both sums o(e''~^)). So, 
from Lemma [3] we get straightforwardly H53p and, from H46p . we get the bound: 



and, when e — >■ 0, we obtain (|54p . 

We may now prove the main proposition of this section: 



□ 



Proposition 1. Consider the stochastic variables defined in i33\} and 128\) and the se- 
quences of positive real numbers {an}^i and {e„}^i satisfying conditions jlSJ) and jUA) . 
i.e. lim af^n^" — 0, for some k G (0, i), and lim a'^e''' = 0, for some (^ € {0, i — 4j)- 

n— >oo ' n— >oo 

Then, ^-alm,ost everywhere w.r.t. sequences of initial data Zoo, 



lim [A„]l sup I 

n^°a te[0,T] 



\Ai%j{t,c)-V^'\t,c)\ Im>, 



= 



where {A„} is an increasing sequence of bounded Lebesgue measurable sets such that _B„ C 
A„ (with Bn defined in (^) and lim A„ — R". 

n—^oc 

Proof. We expand the square power in the expectation value and we calculate the value 
of the three resulting terms. 

Because of definition (|5ip . using (140 p . we can write: 






{fc-2) 



{s,Ch,Zj). 



?{fc) 



We substitute this expression in the definition of yl„ ^ j^j, in system ((33]), and we obtain 
(resumming on the Poissonian variables): 



4e 

„2 Z_^ 



i.j — 1 



[A„]iE" 
dsi ds2 I exp 



aC') 



K,i,Mii,c)\I^M>l} 



[O.t]- 



^J■e[An]L{ I dc 
A 



7^y'"~'>(gl,g2,c) 

[An]L 
ij(k-2). 



-1) 



-Ms[A„]t 



/a dcQJf(sl,S2,c)7^e^' ^(si,S2,c) 



/a„^c 



7^: 



ij {fc-2) 



(S1,S2,C) 



Then, using Lemma |6] we have: 



(57) 



[A„ 



A 



d(fe) 



[An] 



At,c) I. 



{M>1} 



-f f 



dsi ds2 e 



-M[p<'=-2)(si,z0+P<''-^'('i2,2j)(l-<i,j)l 



x(l + e' 



l-d^>j(3)-. 



1 



dcQl^{si,S2,c)'Jt, 



ij (k-2) 



(S1,S2,C) 
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where J?£^ = Re (si, S2, zt, Zj) < KtRi (si, S2, Zi, Zj) is such that 

1 —d n 

J2 \\Rf'^\\L^lo,Tmo,T] < KR.ie^ + oia-")). 






We recall that, using Lemma[6]and since q G ^(R''), we have, ^-almost everywhere 
w.r.t. sequences of initial data Zoc and for any D C R"*, the bounds 



„1 — d n 



(^*) ~^^ ll-Rf'l|z,~[o,T]x[o,T] / dcQ'^{s,t,c)<KqQai\\q\\^{£''+o{a„'^)), 



i,j = '^ 






i.j — l •' U(t.Zt) i — i 

Therefore, using the definition of T^.^"* , given in (|51|l , we get the equality 



[A„ 



I 1 2 

-4„,5,m(*.c) I{m>i} 



= R^f\t)+ 









where 

sup Rt\t) < 4A-,er"||g|Ua^(e'^ + 0(0"'')) + Oie^-^at) -^ 0, 
te[o,T] 

since lim ate'' = with C G (0, i - t^). 

n— foo ^ ^" 

Let 2 = (x, v) and 5 = (j/, w). Define the sequence of measurable sets A^j [z] = {a £ H 
36 G An s.t. a = b — x{s)} and the sequence of functions QXb {z, x) = (/as [z]<ln) * Qnix] 
Since g is a radial function, we have the following equality: 

1 ^ r 

1 y g-M[(p"'-^'(.i...)+P<'=-^'(.2,.,)(i-*„)l / dcQ^^(si,S2,c) = 

/,„,^^,,, dzd^K ®M^](2,2)e-''[''<'"''(^-^)+^"=^''(»-^-)lQ;.,(i,a;(si) -y(s2)) 
where the last term is bounded by 



d 



- fdz^l{q„*q,,){x{s,)-yis2))<2^\\q\\^e. 



2 

n J ■" - . . . . - . . , ^ 

In the first term we can use Fubini's theorem, rewriting 



d2:dz[/iJi ® /i°](2, 2) e"'''' ^"^''''q„{a — x{si))\ x 



da 



X [e-''-<''"''(^-^-)g„(« - y(s2))/Aj;2[„(a ^ ^(.2))] 

and, because of (|42p . we can use Lebesgue's dominated convergence theorem to pass to 
the limit into the integral. 
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{A„}^i is an increasing sequence of Lebesgue measurable set and /a„ — >■ 1 in I/°°(R'*). 
From the J^-a.e. weak convergence 0j^^j^qnlJ-n -^ 0^^iSofo and the L°° convergence of 
/a= [j;] 5 we have 

/[0,t]2x(K<'xR<')2 dsids2dzdz[fj!i (gi nl](z,z)x 

^4,dc|l7W(i,c)|', 
so that we obtain finally 



lim —5-7 / dsi ds2dcQ^^e~ 



i,j = lJ[0,t]2xA, 



dc '[/('=' (i,c) 



|2 



We compute now the double product and, because of the definitions of p^ , An and of 
the boundedness of V'''°\ we get : 

[AuIlE" [ii';^,M(*.c)V^('='(t,c)/{M>i}] ^0{aie''-^) + 

X Xv^ dcV'-''\t,c)qn{xj{si) - c)) . 
Using (|47p . ^-a.e. with respect to Zoo, we get: 

[AnliE" [ii';l,M(t,c)y('=)(i,c)/{M>l}] = 



1 " /■* 



-^cp'■>'-■^>{=l,^i) / j^f/C") 



le""" '^^I'^'M dcV"'>it,c)q„{x,{si)-c) + o{e) 
"'a„ 

so that, because of the ,f^-a.e. weak convergence of qnUn and the hypothesis on /o, we 
obtain: 

-2 lim [A„]iE" [ii''-),,(t,c)V-('='(i,c)/{„>i}] =-2 / dc \v'-'\t,c)\\ 

For the last term we have: 

[A„]iE" [ll/(*'(i,c)l^/{M>i}] =(l-e-'^'"^"'^) / dc\v'^''\t,c)\\ 

'- ^ ■'An ' ' 

I 2 

whose limit is L^ dc V'^'^-' (t, c) because of the L°° convergence of /a„ • 

Hence, we have established Proposition [l] D 

In addition to Proposition [U we shall need its corollary 

Corollary 2. Under the same hypothesis as Proposition\l\ for the obstacles life functions 
defined in i30\) and in i33\) and for all cj) € Ck(K''), ^-almost everywhere w.r.t. sequences 
of initial data Zoo ■ 

M 

lim E"[e'*"^ y^/supp^(ci)jf)„,£,M(i,Ci) - f;(i,Ci)|] =0. 
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Proof. Since for all J > (cf. Remark [T|) 



E"[£''"'X]/supp^(Q)|?)n,.,M(i,Ci)-77(i,C,)|l < 



aW 



-TE."[s''-'Mh^pp4c)\V"'>{t, c) - Al'i j,,{t, c)\] + 2Sfi |supp<?i n A„| , 

o 
the proof follows easily bounding the first term on the right-hand side by Cauchy-Schwarz's 
inequality and using the identity 

E"[{e''-'-Mfl,upp<i,{c)] =^i^[A„]L[snpp(|>nAn]L 
and Proposition [T] D 

We shall need moreover a modified form of the proposition, whose proof will be only 
sketched, being essentially the same as the one of Proposition [T] 

Proposition 2. Under the same hypothesis as Proposttion\ll^- almost everywhere w.r.t. 
sequences of initial data Zoc •' 



lim sup ( -!-E" hAl^l 

n->cx) i<h<„ V \Vh\ L 



\7('='i^(T,„,,,ci)M/,,gr,(T,.,) 



0. 



Proof. The proof is obtained along the same line as the proof of Prop. [T] since Tz^.ci 
is independent of the stochastic times Too, is bounded by T and all quantities involved 
depend on Tz^,ci in a simple way, so to allow to get estimates uniform in Zh. 
Using, instead of (|58p . the estimates 

l-d n 

— 2~ X] ll-R?''^'lli~[o,T]x[o,T]Qn(s,t,c) < Krai{e'' +o{a~'^)), 



«,J=1 



— y2 dcQl^{s,t,c) < KQlvhle"^ ^a^(e'^ + o(a,/)) 

which are a consequence of Lemma [6] and of the bounds (|42p . (|43|l . we perform the change 
of variables (|47|) and we get: 



-E" 



\Vh 

tJ-C'd 



aW 



y2 ds dsids2Qmsi,S2,Xh{s))e 



= 0{at{s^ + o{a-^))) + 

-M[p<'°"^'{si,2i) + P<'=-^'(s2,^i)l 



while for the double product, we have 

^E" \MI^^,^r,^r..,))Al'l,,(t,c,)V^''\t,c^)]=Oiai+'e) + 

Vh I- J 



\Vh 



^ , = l-'0 -'0 



{si,Zj)]. 



Of course 
1 

\vh 



■E" 



Af \/^"Vi,ci) Ic,e%(T,.,) 



^fiCa ds\V^''>{s,xh{s))\ +0{e). 



Then, observing that, thanks to the £P-a.e. weak convergence of qnfJ-„, the quantity 

vanishes for all fixed /i G N, uniformly in Zh because of the choice of the initial data (cf. 
Remark [9] in Appendix I6.2[l and since all error estimates in this section are themselves 
uniform in z^, the Proposition is proved. D 
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4.5. Equivalence between (|33[) and (|31[1 : the O, "P-frozen system. The last step 
to complete the proof of our main theorem is to estimate the distance between the par- 
ticle system 1)33^ and the particle system (|3ip . The difficulty, here, originates from the 
complicate dependence of the life functions, defining both systems, on the stochastic con- 
figuration: in particular, this dependence compels us to use 1)24^ for estimating differences 
of life functions for both species, preventing the use of any of its simplified forms, men- 
tioned in Remark [l] 

Since the system (|33|) is equivalent in the large particles number limit to the system H25|) . 
a mean-field system, what should happen for systems p3|) and (|3fp to be equivalent in the 
limit is that (|31|) is somehow stable (in the large particles number limit) w.r.t. deletion (or 
addition) of a finite number of particles of the two species. This is in fact the case, as we 
shall show in the first Lemma in this section. Then, deleting light particles and obstacles 
in a suitable way from system (|31|) . we shall be able to simplify the dependence on the 
stochastic configuration in the estimates and finally to prove the asymptotic equivalence 
between (l33l) and (|3T1). 



Let O and V be finite subsets of N and k the level of approximation. In this section, 
we introduce a sequence of particle systems where the obstacles with labels in O (resp. 
light particles with labels in "P) have no interaction with the light particles (resp. with 
the obstacles). These systems will give us a tool to estimate the required distance. 

For a given configuration of obstacles, cm ~ (ci, . . . , cm) and a particle initial datum 
Zn = {zi, ■ ■ ■ , Zn), fix the sets of integers O C {f , . . . , M}, V C {f , . . . , n} and define, for 
integers M, k>l: 



ln,£,0 V'jZj) — 7„^^^o V'jZj) — 1^ 



i'^'P}, 



(0,0,-P)/, N 


- hit-P}, 


K,e,M (t,Ci)-- 


= httO}, 


(k.O,V),, X 


= 1 

{^j(s)t U B,(cft)«('=-l-O.T') 


c(k,O..V),, X 


^hA^lf,,rit,c,)<r^^}hiiO}, 



(59) 7;,;,M (*'^j) = ^, ,,."„. ,,,._, ^^^ „ ,/o-^p}' 



K„e,M (*>«) = -E Jo'?"(2:»(s)-Cf)7„,e,M '{.S,Z^)ds. 

Remark that yl^ ^ ^j {t, ci) does not depend on Tc/^ , Xj , Vj for any h £ O and any j € V, 
and it does not depend on Ch, h £ O, whenever (. ^ h. Moreover, for (OyV) — (0, 0), the 
system defined above is the fc-th level approximation system defined in (|3ip . Section [3. 21 
i.e. 

4(fe,0,0) _ w{fc) _,{fc,0,0) _ f(fc) j(fe,0,0) _ „{fe) 



4.5.f. Proof of the equivalence between II33\) and iSlp using the OyV -frozen system. Using 
the bound (|4ip . we can use the same strategy as in [NOR] . 

To shorten the notation, we define in this paragraph, for sets O, 0\ C {1, . . . , M } and 
P,PiC{f,...,n}: 

for M = f , . . ., 
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We prove the following lemmas: 

Lemma 7. Consider the stochastic variables defined in I159\) . Then, under the same 
hypothesis as m Proposition\l\ for all integers k > 1, p > 1 and for all OijVi C N s.t. 
< ^Oij^Pi < CO, ^-almost everywhere w.r.t. sequences of initial data Zoo, vie have 
the limits 

(60) 

1 if/'*' 






xI°.f-r^-*^(Q)M*{»---'-^>-UST'(C?i.^i)m,„.c.,Q) 



= 



Ito sup 1 E"[x,7^.^.-7^-*^(c,)<;''=)(Oi,P0m,„,.,,Q)l =0 



/or all ii, . . . ,ip £ O U Oi and Sjg < T, Zjg — [xj^ , Vjg ) . Limits i60\) are valid also in the 
case where T replaces T^-^a- 

Remark 5. This lemma shows that particle systems obtained from (|3ip by deleting a finite 
number of light particles and/or obstacles are equivalent, in the prescribed asymptotics. 

Proof. The proof is obtained by induction. 
For k — 1 and p > 0, we have 

AoT\OuVi)iT,^^,c,,ct) < ^ J2 r^""' dsqr.{xn{s)-Ci). 
and the second limit in (|60l) follows straightforwardly, since 

'rOUOi,'PU'Pi,M / \ ,nM{l),^^ T, w^ -, , 

T-OUOi.PU-Pi.Af /„ >l"n MT, II II rp 

^e,s,^,.,^ (q) — fPillgll^T 

and E"[2,°,^^.%^^i-"(q)1 < ^iCiT\v,X for p > 0. 

Denoting A^ = #{ii, . . . ,ip,t.}, we build a partition of labels by grouping them in the 
following way. 

We start from the label t and we denote as i (i) , . . . , i (i) , with fcj < fcj < . . . < kwl , 
the wi labels among the p+ 1 labels ii, . . . ,ip,£ having common value £; notice that 
ku>i = P + 1, being c^ ^^ = ci the obstacle associated to the light particle label jp+i- 
We then consider if^ = max{is : is ^ (.} and we call i (2) , . . . , J (2) the W2 labels having 

^1 fciu2 

common value i/j, always using the ordering k\ < k^ < ■ ■ ■ < fc^a (i-6-) '^™2 = .A)- We 
build in this way A'^ groups of labels and we denote q the group label such that fcj' = 1 . 
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We can write then 
(61) ^ T. ^ 



l<jl<n 
l<jp<n 



[0,T]P ;,^i 



xl 



ouOi,vuVi,M^x^if={ii,...,ip,i}-i}_ y^ 

jp+lSVi 



'jp'^e 



dsqnixjp^,{s) ~ ce) 



V^^ E E E •■• E n/ '^^V') / ,. 

^ J\.(i)6'Pi 1<:JJ1)<" 1<JJ2)<" 1<J (n)<"9=1-'^" "'9 -^ IO.tj 



E^M^] 1 



h. — l,...,TJJl — 1 /l — l,...,TJJg h.— l,...,l(;7V 



n 



We use then the bound 

/ 






-fc"' -fc 



/i /i 



</ 



so that, thanks to (|42p . we have 



{ci e%(sj 









[A„i 



AT 



n / n 



JV-l 



H'^C', (,) V^c. ,_,Gre(s,„, 



n 

9=1 



E /_/%(.)-,^"K(^)-,K(.)_,)-%(.)_,)^K(.)^^^'^ 



4"'^!' \<«)- 



and, using (gS]) and E:?[M^] < A'(^^[A„])^ we get finally: 



^<Kl#ri\\q\\ooT- 



i\ (iiTCdK2Y {TK^r+\ 



iTT 



A-o 



so that (|60|) is valid for fc = 1. 
For fc = 2, since 

^OUOi.-puPi.M/ X ^nM(2) 






(c,)^;:;,^^'(Oi,Pi)(r.^-^,.„c,)< 



^OuC»i,T'UT'i,M^ 






- y^ / " dsq„{xh{s) - q) + 



/iGPi 



1 



E 



dsg„(a;h(s)-Q) ^ 'f{c„er,(s,.,z,j} 



in the same way (using Lemma [2] and Lemma [3]) we get the second limit in (|60ll and 

dSH < A^D ( #Pi— + #0(a„e' +e )1( — ^^ 1 ^ — ^^ \vj„\ 

(with iCo depending on |lg|loo, T, Ki, K2, Cd)- 
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Let assume then that (|60|) is true for fc — 1, with k > 2. From the definition of 
■^n s M \^^j .c«i Cf) and the bound (|41|l we have, for p > 0: 

(62) 






s) - Ci)X 



and from Lemma [T] for 5 > 0, 

/ao\ 117 n r('i:-2,OuOi,'PUT'i) r(fc-2,C),-p) I /rri n] ^c,c, 

(63) E,™ [I^^^^.M --^n.^.M |(rz,p^^,c„,c„)J <2(5+ 

iE.„ [Al'^'J'-'\Or U {m},Pi U a+i})(T.,^^,,c„,c„)] 

+ ^E^™ [-4ot^o;7pupj({m},{>+i})(r^^.^^^,c„,c™)] . 

When substituting (|63p through (|62l) into the expectation values in formula (|60p . the 
contribution to the expectation value coming from the first term on the righthand side of 
(p3)) is bounded by 2Ki5, where Ki = Ki{T, Cd, M, Ki, K2,p, #{ii, ■ ■ ■ ,ip, £}), while the 
contribution to (|60|) due to the last terms in (|63p vanish asymptotically because of the 
inductive hypothesis. The remaining term coming from H62p is 0{a'l^n~^). 

The last part of the Lemma is proved replacing Tz^ ,ci by T. D 



Lemma 8. Consider the stochastic variables defined m i33\} and i59\) . Under the same 
hypothesis as in Proposition]^ Vfc > 1, WOjV C N such that < #0, #7^ < <x, ^-almost 
everywhere w.r.t. sequences of initial data Zoo, 

(64) lim sup ^E" [X,^^'r (c.)i<5r' " ^I'ImK^.^.c, , q)|] = 0. 

rn-oo !<„<„ \Vu\ 

Proof. We can prove Lemma [8] by induction. 
For fc = 1 and 1 < ii < n we have: 

I'^u I ^ 

and for k = 2 (using Lemma [S)) 

^li I ^ 

so that (|64l) . because of (|18p and (|19|) . is valid (ii'm depends on ||q||oo, T, ji, K2, Cd)- 
Let assume (|64|l is valid for a given fc > 1 and VO, P. 
We shall use the notation "^2 ~ "12 and "^2 = X] • 

j^P l<j<n h^O \<h<M 

From the definitions (1591) and (I33II we obtain 



r^E"[X,^^-f^(c,)|<f/' -i«,,,i(r.„,.„c,)] < (#p)C,TVllg|k^+ 

jfu I 77, 
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Since both |'*~^'(s,a;j,i'j) and "/'''''^•'^•'^\s,Xj ,Vj) are of the form (gO}, using (|4T|l and 
then Lemma [2] and Lemma [S] ,f3^-almost everywhere w.r.t. sequences of initial data Zoo 
we get (we omit the time integral): 



(65) i^E" 



j^-p 



■tO,V,M f„ \ 



?4^.(^)-c.)|p-^'-7''-^'°-^'|(^,^.) 



rCd^ll(7||ooffl„ 



TCa{#0){- V \v,\) e"-' + SK^e^^ + oia-")) 



+ 






-rO,V,M 

— -7— n 9u(xj(s)-q) 2^ J,^gre(s,3,)k;' -'5n,^,Af (r.,,c;,,c,.. 



(c.) 



We use then lemma[T]and we write, for 5 > 0: 



(66) 



[l'7<^-^^-e:M°'"'lm,.c.,e.)]<25+ 



^[{|^. 



fc-2) 4(fc-2,Ou{£},-pu{j}) , -nAf(fc-2) 



A 



n,e,M 



^ 



O.P 



'({n,{j})}(r.,.c.,c,)], 



since ^i'^gj./^^ (r^^.c^jCh) does not depend on r^. 



The contribution to the right-hand side of (|65|l coming from the last term in the last 
expectation value of (|66p vanishes in the limit thanks to Lemma [7] and the contribution 
coming from the first term is bounded by 25{Ka + -ft'taj^e'' + o{a^'')) for each 5 > 0. 

We evaluate now the remaining term, recalling that h ^ £. 

When h^ i, since bc^-^) - A[^';,"^'^''^*''^''^-''^'| iT,^,c^,Ch) is independent of q, we 
have (again omitting the time integral): 



'(c.) 



\Vu\n 



J2<l"i^jis)-ct)Ic^e%(s,^,)\v^'' ^'-^ 



{fc-2) .(k-2,OU{l},VUi]}) 



iT^j,c^,Ch) 



Jlfp 



< ^^4H^ E E"[^^er.(T,.„)gn(x,(,s)- c)] 



n Dii 



Iffp 



X E" tef (cr) b(^-^) - <-if "^*> 



'^^^^>'|(r.^.,c„co] 



By the triangular inequality, we may then bound the time integral in [0, T] of this 
quantity, ^-almost everywhere w.r.t. sequences of initial data Zoo, by the sum of: 



Me [An 



-^^ / dsE"[/{cgre(T,z„)}gn(%(s 



xE"[l5::f (cO l^^'^-^) -it:llj (T.^.cc,)] 



< A'v,/ sup E" 

l<ii<n 



(fc-2) 



t2:*'t..„(ci) i/c^-^) - <-l,^ (r.„.c,,ci) 
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vanishing because of Proposition [2] and 

xE"fe';';"(ci)M 



-^V / dsE"[/{cgraT,z„)}gn(a;j(s)-c)] 

t-2,Ou{<;},T'U{j})| /^ n] 






< Ka sup E" 

l<j<n 



rO,V,M 






Ak-2,ou{e).vu{j}) _ J(fc-2) L,^ V 



vanishing because of the inductive hypothesis (the constants there depend on Cd, Mi T, 
K\ , K2 ) ■ Since, because of condition a^e^ — >■ 0, all terms vanish, the Lemma is proved. D 

Collecting all results in this section we prove 

Proposition 3. Consider the stochastic variables defined in iSl]) and f33|) . Then, under 
the same hypothesis as in Proposition\7\ ^-almost everywhere w.r.t. sequences of initial 
data Zoo, 



(67) 



1 ^^ 

n— ^00 n ^ — ^ LI ' 



J=l 



and 

(68) hm E" [e''-'A//,,pp4c,)J77« ,,,(t,c,) - C',m(*'=^)11 = 0- 

Proof. We start from (|67p . bounding it first as 



j=i 



j=i 



Using lemma [T] and the bound 



we get, V(5 > 0: 

E.. [l'?i':iM-'7'''l(-,c.)] <25 + iE., [l-4^'/'''({M,{j})l(-,c.) 

4--. [{ 



Define 



,(fc) 



1 T.7irT-{h},{i},A// N .nM{k) 



siy{t) = sup ^E"p.\"f;"^'^"(cO-4„^"'({/.},{j})(r.,,,c, 



+ sup ^E"[xf,,,^(a)ii/"='-ii':^,,,r(r.„.,,Q)] 

l<i<n pi I 

+ sup ^E"[xi,':?;f^>'^^(coi<5i'^^>'-ii':Uim.c„Q)]. 



We get, for aU 5 > 0: 



1 " / 1 " \ 

^EE"[|eM-eMi(M.)]<^Ei-i 



CdrAt5+7 sup Sh,h{t) 



te[o,T] 



and (|67|) follows from Lemma [T] Proposition [2] and Lemma [S] 
As for (|68[). we use Lemma [T] to get the bounds 



[i'?i';U-CMi(i-c,)]<25+ 
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iE.,[{!^^,;^«({M,b-})i + iAh<t>'^^"-A 



n , e , M 



S,«l}(*- 



Cf, 



and, for h ^ i, 



(fc-2) _ rfc-2,{«},{j}| 
n,e,A/ "n,E,M I 



(i,c,0 <25 + 



I^- [■ 



-4; 



ni\/{fc-2) 
0,0 



({M,b-})l(i,c,) , 



and we bound the difference of risk functions through (|4ip 

.riM 

"^0,0 
n 

when ce ^ [J 7l{t,Zi), we get then, through (|69p . for all 5 > 



Noticing that the quantity -Agg ({^},{j})(t, c^) does not depend on ct for h, 7^ 



'[e'^-^i 



E" k''-^M/,„pp4c,)|<i,M(i,c,) - Ci.M(*.c^)l I < ^,[snpp4']LK,{5 + ^) 



'(53 



,(fc-2). 



j,(fc-2). 



sup 4'''f-l'(*) +"n^'' ^ sup d'/'^ ^'(s) 
sS[0,T] ' sg[0,T] 



+e 



d-ll 



V sup -j-rl 



[0,T] 



dsiq„{xi{si) - C£)7{cjg7;(s,2)} 



where Kp = ii'p(r, ||5||oo,-ft'i, -^"2, Mi C'd) and Kd = -ft'd(/i, 0, ||1i;1/o||li)- The proof of the 
proposition follows using Lemmas [71 18] and Proposition [2l D 



5. Final proposition: asymptotic equivalence of ([3T|) and (f25l) 

We may now prove our final proposition, which will allow us to establish the vanishing 
Umits of (El) and (133 : 



Proposition 4. For the life functions defined in i30p and in i31p . under the same hy- 
pothesis as in Proposition\7\ and for all 4> £ Ci,(K'' x M.'') and ip G Ck (R''), the following 
limits are valid, ^-almost everywhere w.r.t. sequences of initial data Zoo'- 



(70) 
and 
(71) 



1 ^ 

lim E"[!^^0(r*(2O)(e<'='(t,^O-CM(*-^O)ll =0 



lim E"[e''-^X^J,.pp^(cO|r?i'=^M(i,cO-r^('='(t,cOll =0. 

Proof. By the triangular inequality: 

1 '^ 

1 "■ ll^ll ^ 

E" [li ^ <p{T\z.))[e' - e,Mi(t, ^01 + ^ E \iiiM - es,Mi(^, -0 

1=1 8 = 1 

and the right-hand side term vanishes because of Lemma [5] and Proposition O 
In the same way, 

M 



A(fe) 



E"[£''-^M7e.pp^(c){|<i,^ - C^.mI + IC.M - '?^'"^l}(i,c)] 
and the right-hand side term vanishes because of Corollary |2] and Proposition |3] □ 
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This proposition completes the proof of theorem [T] 



6. Appendix 



6.1. Existence and Uniqueness of solutions for the limit system. We give here 
the theorem of existence and uniqueness of solutions of system ©, which we can state as 
follows: 

Theorem 2. Letfo>0 and ao>0 be s.t. foeL\Ri-W^'°°{Ri)),vfo€L\Rt;W^'°°(Ri))n 
L°°(m.i X R^), v'^fo G L^{Rt;L°°{M.i)) and oq G W^'°°{Wi). Then for each T > there 
exists an unique solution {f,o-) to the initial value problem {^ in the interval [0,r]. 

The proof of theorem ([2]) is the consequence of the following proposition: 

Proposition 5. Consider the space 

W = {F={Fi,F2,F3) : F^>Qi = 1,2,3, 

F-i G L°°([0,T] xK"" xR''),F2,Fs G L°°([0,T] x R'') } 
with norm 

\\F\\w = l|-P'l|lL~{[0,TlxM'ixK<') + l|-f2|ll,~{[0,TlxR'') + 11-^3 H^oo ([o,T] xK'') 

and the map M = (Mi, M2, M3) : W -5> W defined, for fo,ao > 0, as: 
Mi[F]{t,x,v) = /o(i-'ut,u)e-^''l"l^o'^-'F3(-,— (t-»)) 

(72) M2[F]{t,x) = 4, dvfoix - ^;i,^;)e-C•<il"l/o dsF,(s,.-.it-s)) 

M3[F]{t,x) = ao{x)e''^-f'*''''''^^'-''K 



Assume ao G W^'°°{R') and /„ G L\Ri;W^'°°{Ri)), together with vfo G L\Ri;W^'°°{Ri))n 
L°°(R^xR^) andv^foe L^{Rt;L°°{Ri)). 
Then 

• M IS a strict contraction onW whenT<To, whereTo depends on Cd,®, ||i'/o||Loo(nxR<i)! 

lk/o||Ll(K^;Vl'l,°o(H^))) lkllwl.°°(K<i)) 



Let A/[f] — f be the (unique) fixed point of M , then for each t G [0, T], T <To, 
Mt,-,-) G L' 






'''""*; M/''°°(R^)), 



G L\Ri 
G L\Ri 
G W^-"°(R'') 



W'°°(R°))nL°°(R" xR^), 
i°°(R')) 



Remark 6. The map (|72|) is defined so to represent the solutions to the linear problem 
(I25p (cf. the proof of Theorem [2] on next page), and therefore it does not depend on 
Fi (the sources in the linear problem are indeed F2 and F3). Three components are (of 
course) nevertheless needed to get a map having as unique fixed point the solution to the 
nonlinear problem ([2]). 

Proof of Proposition[Si Consider F, G G W. Then 

||Mi [F] — Ml [G] II ^oo ([o^T]xR''xR'') < C'dT||«/o|| !,== (RdxR'') II -P3 — G3 ll^oo QO^rjxRd) 

ll-'^a [F] - Af2 [G] II ^00 ([o,T]xR'') < C'dr||t;/o II ^1 (ud.^oo (Rd)) II F3 - G3 II ^.^o ([o,T]xR<i) 

II Af3 [F] - M3 [G] II ^oo ([o,T]xR'') < ©r||(To|| i<x, ([o^TJxR'') II -P2 - G2 II L°° ([0,T]xR<i) 
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We can therefore write 

(73) |1M[F] - M[G]\\w < AT\\F - G||w, 

where^ is a constant depending on Cd,0 and on the norms ||f/o||i:,oo(iidxR'')) I1''"o||l<=°([o,t1xK'')) 
lk/o||i,i{E'J;L°o(R<J))- Whenever T < To, where To — ^, M is a, contraction on the 
complete space W. Therefore, whenever T < -j, there exists a unique fixed point 
{f,g,o-) = M{f,g,a) (and of course, g = J dvf). 

The properties of the fixed point functions listed in the second part of the thesis of the 
proposition follow trivially from H72p and the assumptions on /q. D 

Proof of theorem\^ The map M maps the point F in the point Af[F], with Af2[F] = 
JdvMi[F] and (MJF] ,M3[F]) solution to the linear problem ^^, with sources /jj^dw/''""^' = 
F2 and a^''"^' = Fs and initial data f{0,x,v) — fo{x,v), a{0,x) — ao{x). The fixed point 
f of AI is therefore s.t. f — (f, J dvf, a), where (/, a) is the solution of ^ for t G [0, T], 
with T < ^. 

Because of the property of the fixed point (/, J dvf, a), the solution is prolongeable for 
any value of T > 0. □ 

Remark 7. Since /i(i,-,-) G L'(<; W"''°°(R^)), Jdvfi{t,-,v) G Cb{R''). 

6.2. Condition for the weak convergence Qn/Xo ^ <5o/o- We state here a simple 
condition on the growth rate of the generic term of the sequence {an} such that the 
product of two weakly convergent measures converges weakly to the product of the two 
limit measures. 

Lemma 9. Let ^ be a probability measure defined as ipp, with one particle probability 
density /o s.t. 

(74) Jdvfo{-,v)e.y{9.''), 
and q G ^(R"^) a non negative function s.t. 

q{x)dx = e > 0. 



/■ 



Take a sequence of positive real numbers {anjJJ^i (on > 0^ such that for some k G (0, |) 



lim ^ = 0. 

n—¥oo n^ 



2' 
d 



Then, given <j) G diR"^ x R"^), &>-a.s., 
(75) -y^a'iq{a„Xh)<f>(xh,Vh)-^Q fo{0,v)(j>{0,v)dv. 



Proof. We observe first that, since q G o5^(R'*), we can use the identity 

e' '^ q{k)dk 



r / ik-x ■ 



{2-kY 



VA 



where q G ^(R'') is the Fourier transform of q. 

Because of g G ^(R'*) and (|74p , in all calculations below we can apply Fubini's theorem. 
We can then write: 
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1 " f 

-'^a'iq{a„Xh)(p{xh,Vh) - Q / fo{0,v)4>{0,v)dv 



(76) 



(2 






e''''°'^(j){xh,Vh) 



e '^ fo{x,v)<l){x,v)dxdv 



dk 



+ 



a„q(anx)fo{x,v)(l){x,v)dxdv — Q / fo{0,v)(l>(0,v)dv 



We consider then, for h = 1, . . ., the sequence of independent stochastic variables 



'^hW = e''''^''(f>(xh,vh) - 



e^'"" fo{x,v)(j>{x,v), 



s.t. E,g2[!/Jj (fc)] — and Eg2[ji/,j p(fc)] < 2||(;/)||^ and for which we have, for any even j G N, 

:^(E^[l.l^(fc)])4 < :^ 

n2 n2 



(77) Eg. 
We define then 

(78) ujt : 
and we get: 



1 " 



k 
, , , dk\q( — ) 









<C7^l|g||Li+^^ 



(79) 



1 f 

- y^a^(7(a„x,i)(^(a;h,Uh) - e / fo{0,v)(j}{0,v)dv 

+ \ a„q{a„x)fo{x,v)(l){x,v)dxdv — Q fo{0,v)(j>(0,v)dv 

lition of (i 



Since, because of the characteristic function in the definition of a>,t 



by Tchebycheff inequality and ([77)1 . for 5 > 0, 



^{K\>S)< 



i^d^t) ^ K a(j-i)K.-l 



< —^a„n 



5 -- 5""' 
(with K depending on ||<^||oo and ||(7||li); we obtain then, for a £ (0, |) and for all 5 > 0, 



J2 0^{\ujt\ > S) < oo, and therefore uj^ 



^a 



0. 



Since lim — g- = and qn(x) — a„q{anx) is such that qn -^ OSq, we get finally 

n— foo ^ 

-y^a^g(a„a;h)(?!)(a;h,Uh) A^' 6 / /o(0, t;)(;i(0, ?;)di;. 



D 



Proposition 6. In the same hypothesis as in Lenima\^ on a full measure set with respect 
to 0^ 

q„{x)fio{x,v) -^ QSo{x)fo{x,v). 
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Proof. We consider first functions in the separable space Co (R xR ), the set of continuous 
functions vanishing at infinity. Let D be a countable dense set in Co(R'' x M'') and consider 
the two sets 

A = {Zoo ■.y(t)e D / dxdvqn{x)fio{x,v)(t){x,v) ^ e / dvfo{0,v)(f>{0,v)} 

and 

B = {Zoo ■■ / dxdvq„{x)fio{x,v) -> 9 / dvfo{0,v)}. 

Jr''xR'^ JrA 

Because of Lemma (9] both sets have full measure, so as their intersection An B. 
For Zoo € An B the sequence {qnfi^}^^i is a sequence of finite positive measures s.t. 
for all </> G Co(R'' x R'*) 

(80) / dxdvq„fio<p~^e dvfo{0,v)<p{0,v) < ex. 

Jr'^xR'' Jr'^ 

Since the convergence in (|80p is valid also for <^ = 1, on the set AnB, weakly in the sense 
of measure, 

qn{x)fin{x,v) -^ ldSo{x)fo{x,v) 
(see e.g. [MA| . p. 90, theorem 6.8), and the proposition is proved. D 

Remark 8. Since /o G ^(R'' x R"*) and J dvfo G ^(R**), the convergence proved in 
Proposition [6] is valid also for \v\-'qni-i.o, j = 1, 2, . . . (i.e. \v\-'qrifio ~^ I'^P'^o/o)- To show it, 
it suffices to rewrite the proof with u]^ replaced with Iw^i^Jf. Under the same hypothesis, 
it is also possible to prove, by induction, that ®^lj^g„/i,i -^ ^^LiSofo, for A/ = 1, 2 . . .. 

Remark's. Whenever, in addition to the hypothesis in Lemma[9l we have /o G //^(RjJ; VK^'°°(R!^)) 
and v^ fo G //^(RjJ; L°°(R!^)), we can obtain, following the same steps, the weak conver- 
gence of {Ta'"q)nlJ-n and \v\{Ta'"q)nfJ.n, whcrc we define {Ta'"g){x) — g{x + vt + a). The 
convergence is uniform in a G R'', as can be easily checked. 

Moreover, a very simplified form of the procedure allows to prove that, given a limit 
density /o G ^(R'' x R'') and two positive integers P,Q, on a full measure set w.r.t. £P, 
\v\^^J.o -^ \y?fQ and iSil^-^lvkl' fin -^ <Sii^i\v\\fo, for j = 0, . . . , P and i = 1, . . . , g. 
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